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■ Abstract 

\ Color superconductivity is a possible phase of high density QCD. We present 

a systematic derivation of the transition temperature, Tc, from the QCD La- 
CS| ' grangian through study of the di-quark proper vertex. With this approach, 

■ we confirm the dependence of Tc on the coupling g, namely Tc ~ fig^^e~'^^^ , 
, previously obtained from the one-gluon exchange approximation in the su- 
I per conducting phase. The diagrammatic approach we employ allows us to 
I examine the perturbative expansion of the vertex and the propagators. We 

find an additional 0(1) contribution to the prefactor of the exponential from 
the one-loop quark self energy and that the other one-loop radiative contri- 
. butions and the two gluon exchange vertex contribution are subleading. 

PACS numbers: 12.38Aw, 12.38.-t, ll.lO.Wx, 11.15.Ex 
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I. INTRODUCTION. 



Color superconductivity is a possible phase of high density QCD, pioneered by Bailin 
and Love and others, who pointed out that the 3 channel of a di-quark interaction is 
attractive through one gluon exchange. Recently, using an effective four-fermion interaction 
in the superconducting phase, much work has been carried out on examining an energetically 
favored condensate which, for Nj = 3, breaks the original SU{3Y x SU{3)^^ x SU{3)'^^ 
symmetry to its diagonal symmetry This color superconductivity mechanism has 

been called color-flavor locking. 

Several attempts have been made to determine the parameters of the superphase from the 
QCD action |]6HT2| . In all cases one arrives at the same dependence of the zero-temperature 



gap energy on the QCD running coupling constant evaluated at the Fermi energy, namely 

A ~ fig-'e-^/^, (1.1) 

where n is the chemical potential and g the running Yang-Mills coupling evaluated at /i. 
This dependence upon the coupling differs significantly from the BCS case, A ~ fie~'^ , 
due to the long range propagation of magnetic gluons. The coefficient k = Svr^/v^ was 
first derived by Son [0 , and subsequently verified in ^ and [p!0| , p!T| . The latter two of these 



generalized this scaling behavior to the transition temperature and in addition derived its 
ratio to the zero temperature gap. All of these results were obtained from superconducting 
gap equations with one gluon exchange. As such, although they contain the correct leading 
order behaviour that determines the dependence of the gap on the coupling, more detailed 
calculations are required to recover all the leading order contributions to the pre-exponential 
factor. 

In contrast with previous works, we approach the transition temperature from the normal 
phase, with T ^ yU. There are several advantages with this approach. First of all, the 
propagators in the normal phase are not subjected to modifications from the long range 
order, the form of which are Ansatz dependent. This ensures that issues of gauge invariance 
and higher order corrections are relatively simple to handle. Secondly, the integral equation 
for the proper vertex function, which determines the pairing instability, is linear while the 
gap equation in the superphase is nonlinear. Thirdly, the hard dense loop contribution to 
the gluon propagator is free from the Meissner effect. Though it has not been taken into 
account so far in the gap equation calculations in the superphase, the Meissner effect is 
nevertheless expected to alter the pre-exponential factor 0. 

Working directly from the QCD Lagrangian provides a natural framework within which 
to examine the perturbative nature of the theory at high density. Indeed, we shall find that 
interactions of second order, 0{g^), make a leading order contribution to the pre-exponential 
factor and that all higher order contributions are subleading. The sum of these effects 
suggests that a derivation of the transition temperature from the normal phase will not 
only provide a rigorous verification of the results obtained from within the superphase, but 
should also allow a clean and exact determination of the pre-exponential factor. Combined 
with the aesthetic benefit of calculating directly from the QCD Lagrangian, this more than 
compensates for the technical complexity of this approach. 

Starting from the SU (N) QCD Lagrangian, the calculation of the transition temperature 
can be cast easily into thermal diagrams with gauge invariance manifest. Following the 
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formulation developed by Gor'kov and Melik-Barkhudarov for the non-relativistic Fermi- 
gas, which allows for a systematic examination of the higher order contributions, we 
calculate the transition temperature to leading order in coupling and obtain 



nksTc = i^g ^cc'e 
where up to a constant of 0(1), c = 



6N 7i£ 

N+1 g 



[1.2) 



1024\/27r^A^^ in agreement with IHIIl, and the 



exponent, — 



6N 



N+1 g 



in agreement with [^^|-|ri[], are determined by the leading order one 



gluon exchange process. The previously unreported factor 

c' = exp[-4(vr' + 4)(iV-l) 
L lo 

~ 0.1766 for = 3, 



:i.3) 



comes from the logarithmic suppression of the quasi-particle weight in the dressed quark 
propagator. For Landau damping, obtained in the hard dense loop approximation, the 
contribution to the prefactor from two gluon exchange diagrams is subleading in g. This, 
however, is not the case for a hypothetical static screening case where the perturbative 
nature is completely spoiled by infrared log-enhancement in higher orders. 



II. CALCULATION OF THE QCD TRANSITION TEMPERATURE. 

We consider an SU{N) color gauge field coupled to A''^ flavors of massless quarks with 
the Lagrangian density 

C = - i,n,{d, - tgA,)^f, (2.1) 

where F^^ = df^A^ - d^A"^ + gf^^^A^Al, A^ = A^^f" with r the SU{N) generator in its 
fundamental representation. Since the Lagrangian (p.l|) is diagonal with respect to both 
flavor and chirality, the corresponding indices will be dropped below. 

We derive the transition temperature by investigating the onset of the pairing instability 
in the proper vertex function corresponding to the scattering of two quarks at non-zero 
temperature and chemical potential. This vertex function with zero total momentum and 
zero total Matsubara energy, r^^'^*'.^^^'g^(?T,', n|p',p), is shown in Fig. |I|, where n and n' label 
the Matsubara energies iun = ^(^ + |) of individual quarks. Each of the superscripts 
c, which denote color, are associated with a leg. The subscripts s, which label the states 
above or below the Dirac either -|- or — . We flnd it convenient for the partial wave 

analysis to associate the Dirac spinors u^p) and v{f>), which satisfy the Dirac equations 
(74P — ■ p)u{'p) = and (74P — ■ p)v{'p) = 0, to the quark-gluon vertex instead of to 
the quark propagator. Therefore, where s = (si, S2) represents the incoming subscripts and 
s' = (ss, S4) represents the outgoing subscripts, suppressing the color indices and momentum- 
energy dependence, we write 

Ts'-s = U'yis3,p')Us{Si, ~p')T^s,af3UaiSi,p)Up{s2, -p) ■ (2.2) 
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The vertex function T^s^ap is given by conventional Feynman rules, and U{s,p) = u{p) or 
v{—p) for s = + or — , respectively. The proper vertex function satisfies a Dyson-Schwinger 
equation as shown in Fig. 0. This integral equation of Fredholm type may be written with 
all indices suppressed as 

T{n,n\p',p) = T{n',n\p',p) + I ■777^K{n' ,m\p' ,q)T{m,n\qJ), (2.3) 

where F represents the two quark irreducible vertex with all indices defined in the same way 
as for F. The kernel has the explicit form 

K%^:s%\%{n'Mp\q) = fj^^,%^;^^^{n\m\f,q)SsAni\q)SsA-^\ - q), (2.4) 

where we write Ss{n\p) for the full quark propagator with momentum p and Matsubara 
energy iun- The zeroth order quark propagator reads 

Ss{n\p) = -. (2.5) 

and the diagrammatic expansion of F to 0{g^) is displayed in Fig. 0. 

The transition temperature may be obtained from (2^) by examining the Fredholm 
determinant, T> = det(l — K), which is a function of T and fi. Tq is then given by the 
highest temperature at which D(T, /i) vanishes. To demonstrate the gauge invariance of 
this formulation, we need to extend the integral equation ( |2.3| ) to include scattering with 
arbitrary total momentum and energy. Denote by K' the kernel for which either the total 
momentum or the total energy or both are nonzero. Then the Fredholm determinant with 
arbitrary total momentum and energy factorizes as 

V = <lei{l- K)<let{l- K'). (2.6) 

On the other hand. In V is given by the sum of bubble diagrams shown in Fig. ^. They are 
in fact manifestly gauge invariant as may be seen following the argument given in . 

For the rest of this paper, we shall work in Coulomb gauge, in which the full gluon 
propagator takes the form 



D.,j{k,u) = D^{k,u;){6,, - -^), (2.7) 



Dii(k,uj) = D^(k,uj), (2.8) 

and 

Dij{k, uj) = Dj^ik, to) = 0. (2.9) 

Since F corresponds to di-quark scattering it can be decomposed into irreducible rep- 
resentations of SU{N) by either symmetrization [representation 6 for SU{3)] or antisym- 
metrization [representation 3 for SU{3)] among the initial and final color indices, i.e. 

Ff,;^(n',n|p',p) = V2 5^^(^M^*)^^Ff,,(n',n|p",p) (2.10) 
+ V2 5''^'''5''^''Tf,/n',n\f,p), 
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where (■ ■ ■) and [■ ■ ■] denote symmetrization and antisymmetrization with weight one, re- 
spectively, r may be decomposed similarly. Since the Fermi surface has a pairing instability 
in the presence of even an arbitrarily weak attractive interaction, as is also the case for BCS 
theory, we need only focus on the attractive antisymmetric channel for non-Abelian theories. 

Both Tf, ,.{n',n\p',p) and Tf, ,.{n',n\p',p) can be expanded in terms of Legendre poly- 
nomials, i.e. 

rf,^,{n',n\p',p) = J2-fi,^^{n',n\p',p)Pi{cose). (2.11) 
I 

Substituting such expansions into (|2.3|), we find another Fredholm equation satisfied by 

7^,_,(n',n|p',p): 

li',si^',n\p',p) = ^l, {n,n\p',p) + -J2 / dgK^, „{n',m\p',q)-fl„ {m,n\q,p), (2.12) 



where the kernel Kl, ^ has the form 



K,,^,[n , n\p,p) = 2^^(21 + 1) — ^ '^^^^\P^ ^ '^^'^^P' ' (^.13) 

We consider the / = term in the partial wave expansion, although the higher partial wave 

terms may contribute 0]. 

The Fredholm determinant of ( |2.12D with I = can be written as "D = Y\j{^ — Aj^) with 
the eigenvalues (labeled by the integer suffix j) defined by the homogeneous equation, 

f,,{n,p)=X'-Y^J dqKl^XnMP.Q)fs{m,q). (2.14) 



At sufficiently high temperature, all A| > 1, so that P 7^ and there is no instability — 
the theory is in the normal phase. As the temperature is reduced, we find the transition 
temperature to the superconducting phase is that at which the smallest of {A|} reaches one. 
The solution of (p.l4|) provides the eigenvalues in terms of the parameters of the theory; the 
temperature, coupling and chemical potential. Hence the inversion of Xl{T,g,n) = 1, where 
Ao is the smallest eigenvalue, yields the transition temperature Tc. 

In the presence of a Fermi sea, hard dense loops have to be included in the gluon prop- 
agator. As a result, the Coulomb interaction is strongly screened by the Debye length, 
Xd = 'm~D, where 

TT^ 70 ^^^e/^(''-^) + l - 27r2 • ^^-^^^ 
The dressed Coulomb propagator at momentum-energy (fc, ioj) reads 
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with a^{k,uj) ~ mjj for lu <^ k <^ fi. However, the magnetic interaction is poorly screened. 
While a magnetic mass may exist, of order T, Landau damping |15| prevails at /i ^ ksT. 
In this case the propagator for a magnetic gluon is 



D'%k,uj) 



fc2 + cj2 + a*^(A;,cu) 



(2.17) 



The only region of infrared sensitivity on the (fc, co')-plane is <^ /c -C where a^^{k^uj) ~ 

n', n|p', p) arises from the one-gluon exchange 



|i| ■ '^'^ order, the contribution to 
diagram of Fig. ^ and is given by 

2 . 1 
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12p'p 



In 



\P' -P\^ + f^^D^n' 



3, V 



(2.18) 



where a term finite in the limit g —>■ has been dropped. To the leading order of In fi/kBT, 
the summation over Matsubara energy can be replaced by an integral and \p' — ignored. 

Furthermore, the components /+_, /_+ and / can be neglected and the integral over q can 

be carried out for a solution smooth in the neighborhood of the Fermi sea, p c:^ fi. Equation 
( p.l4| ) is then approximated by 



where 



v' 



5/2 „5 



V 



and 



1024y27rV 



N'/^g^kBT 



(2.19) 



(2.20) 



1024v^7rV' 



(2.21) 



The reduced kernel K,{v^v') is given by 



\' g' 

2 247r2 



\U — W 



+ In 



+ i>'\ 



(2.22) 



with an ultraviolet cutoff z/q ~ 1 introduced. The eigenvalue problem ( |2.19|) can be solved 
using the same approximation employed by Son 0, which amounts to replacing the kernel 
of ( p.l9| ) with (2/z/') lnl/z/> where z/> = max(i>, 0'). We find that 



x2 9 

' 247r2 



1 



In^ 



TT 



where j is an integer and 



Inl/e 



sm 



J + 



1\ lnl/z>' 

TT- 7— 

Inl/e 



(2.23) 



(2.24) 
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Setting the smallest eigenvalue to one, namely Ag = 1, we finally arrive at the exponent and 
the prefactor c of ( |1.2|) . 

Here we wish to highlight the mathematical structure of ( |2.3| ) which characterizes the 
long range attractive interaction. If, instead, the pairing interaction was of a short range 
nature, the transition temperature could be located by means of the standard expansion of 
the Fredholm determinant 



dp K{n,p\n,p) 



+ 



1 



E 

n,n' 



dp / dq 



K{n,p\n,p) 
K{n', q\n,p) 



K{n,p\n', q) 
K{n', q\n', q) 



+ ■ ■ 



(2.25) 



The m*^ term in the expansion on the right hand side contains m kernels folded together and 
for a short range interaction would be of the order g"^"^ In - since there is only one eigenvalue 
of K which diverges as In ^ in the limit e — > 0. At the transition temperature one would 
have g'^ln- ~ 1, and thus the order of magnitude of the subsequent terms would be ^2(m-i)_ 
Therefore only the first two non-trivial terms (m = 1,2) would be sufficient to determine the 
transition temperature up to the leading order of the pre-exponential factor. On the other 
hand, for the present long range attraction of the QCD model, the logarithm in the kernel 
of ( |2.14D makes the m^^ term of the expansion ( p. 251) of the order g"^"^ In^*" - since there are 



now an infinite number of eigenvalues of K which diverge as In - in the limit e — 0, as 



indicated in (p.23|) . Hence the series can not be truncated at Tc, and a new method has to 
be employed to fix both the exponent and the prefactor. 

Note that the second term on the right hand side of ( p.25| ) is in fact logarithmically 
infinite for the kernel ( [2.1^ ). On the other hand, this term corresponds to the sum J2^j'^ 
which is convergent according to ( p.23[ ). The reason for this apparent paradox lies in Son's 
approximation, which led to the eigenvalues given in ( p.23|) . When the approximation is 
corrected, the eigenvalues instead become 



-2 



247r4 



1 + 



In^ 



(2.26) 



with Cj ~ j for j ^ 1 ||T6[; which explains the appearance of the logarithmic divergence 



in this term when summing over j. 

It is interesting to note that a similar non-BCS scaling behavior of the type indicated 
in ( |1.1|) was obtained via a mean field calculation of 2D superconductivity at the von Hove 
singularity of the electronic density of states . 



III. HIGHER ORDER CORRECTIONS. 



We now come to the question of higher order corrections to the kernel from the per- 
turbative expansion of the quark propagator and the irreducible vertex F. These become 
important if sufficient powers of the infrared logarithm accompany the coupling constant g^. 
It follows from ( |1.2| ) that In n/ksTc ^ 1/g at the transition temperature. Thus, if the 0{g'^) 
contribution to Kg, g is of the form g^ln^ n/lun' — Vn\ or g^hi^ ix/\vn'\t its magnitude relative 



7 



to the 0{g^) term given by (|2.18|) will be g^~^. The perturbative robustness of the exponent 
then requires that 6 < 3; robustness of the prefactor requires 6 < 2. We find with Landau 
damping that 6 = 2; one logarithm originates from the leading order of 7, ( p.l8| ), and the 
second from the self-energy of quarks 0. Therefore there will be an 0(1) correction to the 
prefactor. 

Parametrizing the dressed quark propagator above the Dirac 



we find, to one loop order, that 



d 



— ReE(po,p) 
dpo 



g^ 



In 



N 247r2 m'jjmax{\po\,kBT)' 



(3.1) 



(3.2) 



while ^ReE(po,R) remains finite in the limit po ^ 00. For a Matsubara energy po 



i^(n + i), we thus have 



2^' 



S+{n\p) 



l + ^£^\nj^^]u.-p 



(3.3) 



/i 



Following the steps which lead from ( |2.16| ) to ( |2.22| ) above, the inclusion of this effect 
amounts to replacing the reduced kernel /C(z/, u') by /C(z/, u') + A/C(i/, u') where 



A/C(z/, u') 



N 247r2 



/C(z/, z/ ) m — . 



(3.4) 



Treating A/C(i^, z/') as a perturbation, the shift of the eigenvalue in ( p.23| ) with j = turns 
out to be 



5 



1 

^0 



^0 dv r'^o dv' 

^ ^Uv)AJC{v,v')Uv') 
v' 



V 

2{tx^ 



4) 



l\ -I 



9' 



N 



247r2 



In'-. 

e 



(3.5) 



The condition for the critical temperature, Aq = 1, now becomes 



N 



\^247r2 / 



N 



N 



(vr^ + 4)ln^ - = 1, 



(3.6) 



the solution of which, for small g, gives rise to the result ( |1.2D with both c and c'. The 
logarithmic dependence of (|3.2|) upon the coupling constant g, written in m'jj, will change 
the prefactor g~^ of (T72) to g~^+^(3)_ However this correction is of higher order. 

Other higher order corrections to 7 have also been partially addressed in the literature. 
The vertex correction has been discussed in and some renormalization group arguments 
have been applied to the straight box diagrams of two gluon exchange, [0, which was conjec- 
tured to be of the same order as the crossed box diagram. For the case of Landau damping, 
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our analysis of the vertex correction is in agreement with indicating an 0{g) contribution 
to the prefactor. We also find that the crossed box diagram is free from any logarithmic 
enhancement and so its contribution is suppressed relative to the one gluon exchange by 
a factor of g"^/ \n{n/kBT) On the other hand, the straight box diagram, which corre- 
sponds to the convolution of two single gluon exchange kernels, is logarithmically divergent 
at T = and contains all the powers of logarithms at T 7^ necessary to produce the result 
of the ladder sum implicit in ( |2.14| ) . 

The crossed box contributes to 'y^_^_ a term 

1 4 /"i r°° duj r d^l ,r Q rr Q 

B = --9 |/cos^y_ -y j^,D,Ai - -^MD.Ai + -^M 

X [m(p')7m'S'f(P + /, iuj)'-fyu{'p)][u{-p')'-iy: Sf{-P + /, iuj)'-f^,u{-p)], (3.7) 

where |p| = = /x, P = \{p + 'p'), q = p — p' and the summation over Matsubara 
energy has been replaced by an integral for T <^ ^. Ignoring the Coulomb propagator, the 
contribution from the small scattering angle, 16*1 < ^ 1) and infrared region, IR : |Z | -C 
and \uj\ <^ fi, 

^ 2^7-00 yiR27r(27r)3 ^ 2' ^ ^ 2' ^ 

x[u{p')-fiSF{P + l',iujhMp)M-p'hfSF{-P + tiujhi'u{-p)], (3.8) 

is bounded: Pir < b where 

' d9 / dpd¥-^ , , ' , ,^1,, , , (3.9) 



327rV^o JiR {rl + i^\p\){ri + k\p\){p^ + El){p'^ + El) 

Here E± = \P ± I \/ p — 1, r± = |/ ± q\/p, p = \^\/p and k = j-j^- Transforming the 
integration variables from 6, r to E±, r±, we end up with b = ^^^4^2 Jo dpK{p) where 



K{p) = I dE+dE_drldrl J ,^3 , T.2^, ... , z.2 ^ ^ (3-10) 



rj^r^\E^E^ — p^' 
{r% + K\p\){rl + K\p\){p^ + El){p-' + Ely 

with the Jacobian 

J = [(E+ - E_Y - A{rl + )(E+ - E_f - 16(E+ + E_f + IGr^ ]-V2. (3.11) 

As p 0, we find that K{p) — > const xp^^/'^ up to some power of In p. Therefore Bm as 
well as B is free from infrared divergences. 

For the sake of comparison, we have also examined the 0{g^) corrections with only a 
static mass gap for gluons, m -C /x. This amounts to replacing the magnetic gluon propagator 
(Qby 

D'\Kuj) = -^ (3.12) 
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The integration over the Euchdean energy can be carried out readily by the residue theorem 
and the remaining integral can be classified according to the contributions from the gluon 
and quark poles. Denoting the 0{g'^) contribution to by A7, we find to the leading 
order in In/x/m 



A7 



27r2 



2N\ N 



1 ^Ar 2 1 



C2 



In^ii, 

m 



(3.13) 



where Ci corresponds to the vertex corrections and C2 corresponds to the crossed box diagram 
in Fig. ^ We write ci = gi + qi and C2 = (72 + where gi and g2 come from the gluon poles 
and qi and q2 come from the quark poles, respectively. The g and q coefficients are tabulated 
in Table H for various cases. This presence of fi/m will ruin the perturbative nature of the 
formulation, provided m T. At this point, the difference between Landau damping and 
static screening is clear. To the leading order of one gluon exchange, the difference merely 
amounts to the substitution of the static screening mass by m^^|a;|^/^, with uj the Euclidean 
energy transfer and |ci;| <C /x. This is not at all the case for higher order corrections, including 
two gluon exchange, even though the infra-red sensitive region for the loop momentum of 
the quark pole contribution in the case of static screening coincides with that of Landau 
damping. 



IV. CONCLUSION. 

In conclusion, we have derived the superconducting transition temperature with thermal 
diagrams in the normal phase. This ensures gauge invariance to all orders. We have also 
examined systematically the 0{g^) corrections and found an additional contribution to the 
pre-exponential factor of Tc in the literature. Unlike the situations with static screening, 
Landau damping significantly improves the infra-red behavior of the higher order diagrams 
and makes the perturbative expansion of the prefactor in terms of g legitimate. 



In a forthcoming publication |16|, we develop a perturbative formulation that enables us 
to eliminate the 0(1) uncertainty of the prefactor c of Eqn. ( p. . 1|) . Since the perturbative 
expansion is given in terms oi g 01 g In g, the application to realistic situations is in fact rather 
limited. As was pointed out in 0, even with asymptotic freedom, g = 0.67 at /i = 10^°Mev, 
and yet the energy scale probed at RHIC is only a few hundred Mev. Nevertheless, it reveals 
some novel properties of superconductivity induced by a long range interaction which has 
not yet been fully examined in the literature. 
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TABLES 



transverse 



covariant 



91 
Qi 

92 
92 



with FS 





1/16 



without FS 

-1/4 






with FS 




-2/3 
1 



without FS 



-2/3 




TABLE I. Leading log coefficients 
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FIGURES 




FIG. 2. The Schwinger-Dyson equation. As in Fig. |T|, F is represented by singly hashed vertices 
and F is represented by double hashed vertices. The expansion of F is given up to 0{g^). 




FIG. 3. The Bubble Diagrams. 
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